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Ultracold collisions between spin-polarized Na atoms and 
vibrationally excited Na2 molecules are investigated theoret- 
ically, using both an inelastic formalism (neglecting atom ex- 
change channels) and a reactive formalism (including atom 
exchange). Calculations are carried out on both pairwise 
additive and non-additive potential energy surfaces for the 
quartet electronic state. In both inelastic and reactive cal- 
culations, the Wigner threshold laws are followed for ener- 
gies below 10~ 6 K. It is found that vibrational relaxation 
processes dominate elastic processes for temperatures below 
10~ 3 — 10 -4 K. For temperatures below 10 -5 K, the rate co- 
efficients for vibrational relaxation (v = 1 — ► 0) from the full 
calculation are 4.8 x 10' 11 and 5.2 x 10~ 10 cm 3 s _1 for the 
additive and non-additive potentials respectively. 



Methods for creating diatomic molecules in atomic 
Bose-Einstein condensates (BECs) are now starting to be 
realised experimentally. They include photoassociation 
spectroscopy |i"|-fH and magnetic tuning through a Fes- 
hbach resonance . An important long-range goal of 
such experiments is the production of a stable molecular 
BEC. However, in most cases the molecules are produced 
initially in vibrationally excited states, and their lifetime 
is limited by collisions with other atoms or molecules. 
Since the magnetic trap depth is typically 1 mK, any vi- 
brationally or rotationally inelastic collision will release 
enough kinetic energy for both collision partners to be 
ejected from the trap. 

Very little is known about vibrational relaxation rates 
for ultracold alkali dimers. Wynar et al. || produced 
87 Rb2 molecules in very high vibrational levels of the 
ground electronic state in an atomic BEC by stimu- 
lated Raman scattering. They analyzed their line shapes 
to obtain an upper bound on the inelastic rate coeffi- 
cient fc lncl < 8 x 10 -11 cm 3 /s. In subsequent experi- 
ments |5| on ultracold (but not condensed) Rb2 molecules 
in a different set of vibrational levels, they measured 
k mcl = 3 x 10~ n cm 3 /s. These rates are too high for the 
production of long-lived molecules, but there is a hope 
that lower-lying states will relax more slowly ||[Io| and 
that methods for stabilizing the molecules can be found. 

Several calculations have been carried out on the vi- 
brational relaxation of molecules such as H2 IqilQll and 



CO Jll|,^2| at ultralow energies. Spin-changing collisions 
of O2 have also been investigated. p3| , [[4[ However, colli- 
sions of alkali dimers present new theoretical challenges 
that are not present for collisions of stabler molecules. 
In particular, the potential energy surfaces are such that 
barrierless atom-exchange reactions can occur; even if the 
products are indistinguishable from the reactants, the re- 
active channels must be taken into account in a full treat- 
ment of the collision dynamics. Barrierless reactions are 
significantly different from reactions such as F + H2 , (Tif l 
which have substantial barriers, and have not yet been 
investigated at ultralow energies. This Letter presents an 
initial investigation of alkali + alkali dimer collisions, for 
the case of Na + Na2 collisions occurring on the lowest 
quartet surface for Na3. This corresponds physically to 
collisions of atoms in their "stretched" spin states, with 
F = F max = I + S and \M F \=F. 

In this work we use the potential energy surfaces of 
Higgins et al. [0 for the 1 A A' 2 state of Na 3 . Their full 
(non-additive) potential was obtained from a grid of ah 
initio calculations with a large basis set, interpolated us- 
ing the reproducing-kernel Hilbert space scheme of Ho 
and Rabitz The potential has a global minimum 

at —1222.1 K, with the atoms in an equilateral triangle 
configuration 4.41 A apart. The symmetric linear ge- 
ometry is a saddle point at —554.4 K (667.7 K above 
the minimum), with r = 5.10 A. Comparison with the 
corresponding triplet Na2 pair potential pj| ], which has 
r e = 5.192 A and D e = 255.7 K, shows that large non- 
additive effects are present. In the present work, we test 
the sensitivity of our results to the potential surface by 
performing calculations using both the full non-additive 
potential and a potential obtained by neglecting the non- 
additive part. 

We initially carried out inelastic scattering calculations 
of vibrational relaxation from the v = 1 state of triplet 
Na 2 , using the MOLSCAT program Such calcula- 

tions are carried out in Jacobi coordinates (R, r, 9) for a 
single arrangement of the atoms, and neglect the reactive 
channels. The wavefunction in the interaction region is 
expanded in a basis set of Na2 vibrational functions that 
are eigenfunctions of the Hamiltonian for free Na2, sup- 
plemented by a wall at large R to provide a representation 
of the Na2 continuum. 

For full close-coupling or coupled states calculations, 
it proved impossible to converge the basis set of Na2 
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vibrational functions. This may readily be understood 
in terms of the potential energy surface. Consider Na 
approaching Na2 for a T-shaped geometry, 9 = 90°, al- 
lowing the Na2 bond length r to relax to minimize the 
energy at each intermolecular distance R. The system 
passes through the minimum-energy equilateral geome- 
try, and then the two atoms of the Na2 move apart to 
allow the Na atom to insert between them. However, 
even at R — (the linear geometry), the energy is be- 
low that for separated Na + Na2. At this configuration, 
the optimum value of r is 10.4 A, and the basis set of 
Na2 vibrational functions centered around r e — 5.1 A is 
completely inadequate for representing such dramatically 
expanded geometries. 



cross sections for collisions of Na2 initially in v = 1 are 
shown in Figs, [l] and || for both the pairwise- additive and 
non-additive potentials. 

The cross sections all follow the expected Wigner 
threshold laws ||,|o]|2l) at energies below 10 6 K. For 
clastic collisions, the cross sections are proportional to 
E°, E 2 and E 3 for s, p and d waves respectively. The E 3 
dependence for the d wave arises because the long-range 
dispersion (R~ e ) term in the atom-diatom potential mod- 
ifies the threshold behavior for I > 1 [^lj . For vibrational 
relaxation, the cross sections are proportional to E~ 1 / 2 , 
E 1 / 2 and E z / 2 for s, p and d waves respectively. The 
oscillations that occur for E > 1 mK arise simply from 
zeroes in sin 2 77, where rj is the scattering phase shift. 




10"" 10 



10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 
Collision energy (K) 



10"' 10" 




10" 



Collision energy (K) 



FIG. 1. Contributions to elastic cross sections from s, p 
and d waves for collinear Na + Na2 (v = 1) collisions, neglect- 
ing reactive channels, (a) additive potential, (b) non-additive 
potential. 

This problem does not arise for collinear scattering, for 
which the vibrational basis set is adequately converged 
for u max = 20. We therefore carried out collinear cal- 
culations of the elastic and vibrational relaxation cross 
sections. The coupled equations were integrated from 
R = 2.3 A to 2500 A. The resulting elastic and inelastic 
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FIG. 2. Contributions to inelastic cross sections (to pro- 
duce Na2 (y = 1) from s, p and d waves for collinear Na + Na2 
(v = 1) collisions, neglecting reactive channels, (a) additive 
potential, (b) non-additive potential. 

The elastic cross sections are dramatically different for 
the additive and non-additive potentials (by a factor of 
almost 10 3 ), while the inelastic cross sections are much 
more similar (differing by a factor of only 10). The large 
difference in the elastic cross sections presumably occurs 
because the zero-energy elastic cross section may be writ- 
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ten as rjg 1 = 47r|a| 2 , where a is the atom-diatom complex 
scattering length |J. This is accidentally close to zero 
for the non-additive potential: a = (—0.90 — 1.19*) nm, 
compared to (33.83 — 0.46z) nm for the additive poten- 
tial. This has the effect that a cl > cr incl for T > 10~ 10 K 
on the additive surface, but only for T > 10~ 4 K on the 
non-additive surface. 

An intriguing feature of the results is the appearance 
of a second linear region (in the log-log plots) for inelas- 
tic s-wave scattering between about 10 -5 and 10 -2 K. 
The cross section in this region is proportional to E^ 1 
for the non-additive potential but to about E~ 4 / 3 for the 
additive potential. A possible reason for this is that the 
vibrational coupling at long range decays as R~ 6 for the 
full potential but as R~ 8 for the pairwise-additive poten- 
tial. 

As described above, the MOLSCAT calculations are re- 
stricted to collinear geometries and a single arrangement 
channel. To go beyond this requires a reactive scatter- 
ing formalism in which all three arrangement channels 
are included. Such calculations are more expensive, and 
have not previously been attempted for systems contain- 
ing three non- hydrogen atoms. 

As a first step towards a rigorous treatment, we have 
performed three-dimensional quantum reactive scatter- 
ing calculations for total nuclear orbital angular momen- 
tum J = 0. The configuration space is divided into an 
inner and an outer region, depending on the atom-diatom 
distance. In the inner region, we use a formalism based on 
body-frame democratic hyperspherical coordinates p3| . 
This has already proved successful in describing atom- 
diatom insertion reactions such as N( 2 D) + H 2 ^ NH 
and 0( 1 D) + H 2 -> OH + H (25pl. A re- 
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lated approach has also been used in studies of three- 
body recombination of ultracold atoms 27 1. The scat- 
tering wave function is expanded on a set of hyperspher- 
ical adiabatic states. This yields a set of close-coupling 
equations, which in our method are solved using the 
Johnson-Manolopoulos log-derivative propagator C^j . In 
the outer region, we use the Arthurs-Dalgarno formalism 
p9| , which is based on Jacobi coordinates. Matching of 
the wavefunctions in the inner and outer regions is per- 
formed on a boundary which is an hypersphere of radius 
25 A. This yields the reactance ^-matrix and the scat- 
tering S'-matrix. 

The inner region starts at a hyperradius of 4 A and 
is split into 297 sectors. The adiabatic states in each 
sector are obtained by a variational expansion on a ba- 
sis of hyperspherical harmonics with A\ symmetry. They 
arc fully symmetric with respect to particle permutations 
to account for the indistinguishability of atoms. At large 
hyperradius, the adiabatic states concentrate into the ar- 
rangement channels and describe Na 2 molecules in even 
j states. At small hyperradius, they span a large fraction 
of configuration space and allow for atom exchange. The 
hyperspherical harmonic basis is truncated at A max , the 



maximum value of the grand angular momentum. A max 
varies from 198 (867 harmonics) at small hyperradius to 
398 (3400 harmonics) at large hyperradius. A fixed num- 
ber of 135 adiabatic states is used in the close-coupling 
expansion in each sector. At the boundary between the 
inner and outer regions, the adiabatic states are projected 
onto a set of Na 2 rovibrational states, with j max = 48, 44, 
40, 36, 30, 26, 20, 10 for vibrational levels v = 0, 1, . . . , 7. 
The boundary between the inner and outer regions was 
placed at a distance such that couplings due to the atom- 
diatom residual interaction can be neglected outside the 
boundary. In the outer region, regular and irregular solu- 
tions of a radial Schrodinger equation which includes the 
isotropic (R~ 6 ) part of the interaction were integrated 
inwards from very large distances (5000 A). 
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FIG. 3. Cross sections from quantum reactive scattering 
calculations for Na + Na2(« = l,j = 0) (s-wave scattering). 
Elastic and quenching results are shown as solid and dotted 
lines, (a) additive potential, (b) non-additive potential. 

Fig. U shows the cross sections as a function of collision 
energy for the additive and non-additive potentials. All 
collisions that are energetically elastic, with or without 
atom exchange, are included in the elastic cross section. 
All other processes (which produce Na 2 (v = 0, j)) con- 
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tribute to the quenching cross section. Na2 rotational 
levels up to j — 20 are energetically accessible at the en- 
ergy of the v = 1 state (23.5 cm -1 ), and all accessible 
levels are populated in the products. 

It may be seen that the Wigner threshold laws for the 
elastic and quenching cross sections hold below 10~ 5 K 
for both potentials. The cross sections are larger for the 
non-additive than for the additive potential, by a factor 
of about 10 for both elastic and quenching cross sections. 
The E~ x / 2 dependence of the quenching cross sections 
corresponds to a constant rate coefficient below 10~ 5 K, 
which is k mcl = 4.8 x 10~ n cm 3 s _1 for the additive 
potential and 5.2 x 10~ 10 cm 3 s" 1 for the non-additive 
potential. The corresponding scattering lengths are a = 
(2.26 - 0.80i) and (2.38 - 8.54«) nm, respectively. 

Outside the Wigner region, the cross sections have 
a more complicated energy dependence. The quench- 
ing probability increases with increasing energy and ap- 
proaches unity at the limit of the Wigner region. The 
quenching cross sections thus vary approximately as E~ x 
above 10~ 4 K, because of the k~ 2 factor in the expres- 
sion for the cross section. Above 10 -3 K, the elastic 
cross sections show oscillations which are similar to but 
less pronounced than the ones in Fig. [TJ. 

Finally, the ratio of quenching to elastic cross sections 
is larger than 1 at energies below 10 -4 K for the additive 
potential and below 10~ 3 K for the non-additive poten- 
tial. It increases up to 500 in the nK range for both 
potentials. 

In future work, we intend to investigate the depen- 
dence of quenching rates on the initial vibrational quan- 
tum number and to investigate the effects of magnetic 
fields and nuclear spin coupling. 

The three-dimensional quantum dynamical calcula- 
tions were performed on a NEC-SX5 vector super- 
computer, through a grant from the "Institut du 
Developpement des Ressources en Informatique Scien- 
tifique" (IDRIS) in Orsay (France). PS is grateful to 
the EPSRC for a Research Associateship under grant no. 
GR/R17522/01. 
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